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Abstract 

Given a Hamiltonian with a continuous symmetry one can generally factorize that sym- 
metry and consider the dynamics on invariant Hilbert Spaces. In Statistical Mechanics this 
procedure is known as the vertex-IRF map, and in certain cases, like rotational invariant 
Hamiltonians, can be implemented via group theoretical techniques. Using this map we 
translate the DMRG method, which applies to Id vertex Hamiltonians, into a formulation 
adequate to study IRF Hamiltonians. The advantage of the IRF formulation of the DMRG 
method ( we name it IRF-DMRG), is that the dimensions of the Hilbert Spaces involved in 
numerical computations are smaller than in the vertex-DMRG, since the degeneracy due 
to the symmetry has been eliminated. The IRF-DMRG admits a natural and geometric 
formulation in terms of the paths or string algebras used in Exactly Integrable Systems 
and Conformal Field Theory. We illustrate the IRF-DMRG method with the study of the 
SOS model which corresponds to the spin 1/2 Heisenberg chain and the RSOS models with 
Coxeter diagram of type A, which correspond to the quantum group invariant XXZ chain. 



PACS numbers: 05.50.+q, ll.lO.Gh, 75.10.Jm 







Introduction 



The Density Matrix Renormalization Group is a powerful numerical real space RG method 
introduced by White in 1992 to study quantum lattice Hamiltonians of interest in Condensed 
Matter and Statistical Physics |l]]. This method has its roots in Wilson's solution of the Kondo 
problem [Q], but it is not confined to impurity problems. The DMRG method overcomes the 
problems of the old Block RG method of the SLAG and Paris groups [^], which in many 
cases gives qualitative correct results but lacks numerical accuracy ( for a review of the Block 
RG method see P, 3) . The DMRG is well suited for Id problems as spin chains 0, |] , but it 
has also been applied successfully to ladder systems and large 2d blocks [|10] . Another related 



developments are : Generalization of the DMRG to classical systems and its relation with the 



Baxter's Gorner Transfer Matrix [11 1, Variational Formulation of the DMRG ground state wave 



function [|T^, Momentum space DMRG |T^, application of the DMRG to transfer matrices |T^, 
DMRG study of quantum systems at finite temperature [T^ , analytic formulation of the DMRG 



IB| . The correlation between blocks inherent to the DMRG method has also been implemented 



in the old Block RG method in references [0, 0, |T^, etc 



The purpose of this paper is to generalize the DMRG to a class of models commonly known 



in Statistical Mechanics as Interaction Round a Face (IRF) or more simply Face models 
In these models the lattice variables are labelled by the points (heights) of a graph Q, and such 
that heights located at nearby sites on the lattice must also be nearest neighbours on the graph. 
The SOS models and the RSOS models, which are a restricted class on the former ones are 
the most interesting examples of IRF models, due to their connection with Integrable Models 



2l| , Affine Lie algebras |2^, Towers of Multi-Matrix Algebras p3|, Q and Gonformal Field 



Theories (GFT) Other important class of models is given by the vertex models, where the 
nearby lattice variables are independent, although the Bolztmann weights may satisfy certain 
conservation laws The well known Heisenberg, t-J and Hubbard models, are Hamiltonian 



or transfer matrix versions of vertex models, which can be studied using the DMRG. It is for 
the later class of theories that applies the standard DMRG. We shall propose in this paper to 
translate the "vertex language", which is used to formulate the DMRG, into "IRF language". 
This translation process is suggested by the fact that some models, like the Baxter's 8 vertex. 



can be mapped into IRF models ^7\, |2^, ^yj. Moreover, if the vertex Hamiltonian has a sym- 
metry described by some group (or quantum group), then the vertex-IRF map consists in the 
factorization of that symmetry. In symbolical terms we may write. 

Symmetry ^ ^ 

In the case of the Heisenberg model, the factorization of the rotational symmetry has lead 
us to formulate the DMRG in IRF variables. The vertex-IRF map is given in this example by 
the tensor product decomposition of irreps of SU{2). The heights coincide with the irreps of 
this group. From (||) it is clear that an advantage in working with IRF variables is that the 
symmetry present in the vertex Hamiltonian is factorized, and consequently the dimension of the 
Hilbert spaces involved is much lower. For numerical purposes this property is also important 
since it implies a reduction of the computational complexity of the problem. On the other hand 
the IRF formulation of the DMRG is the most natural one to discuss its relation with the corner 
transfer matrix formalism and GFT. 

The organization of the paper is as follows. In section I we review the basic concepts and 
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tools of the IRF models. In section II we introduce the real space renormalization group method 
applied to IRF Hamiltonians. In section III we define the density matrix for IRF states and 
use it to propose the IRF-DMRG algorithm. In section IV we apply the IRF-DMRG to the 
SOS model, which corresponds to the spin 1/2 Heisenberg model, and to the RSOS models 
whose graph is a Coxeter diagram of type A, and which are related to the minimal models of 
Belavin, Polyakov and Zamolodchikov [3^. In section V we describe the vertex-IRF map for 



Hamiltonians with spin rotation symmetry, and derive the IRF-DMRG from the vertex-DMRG. 
Finally we state our conclusions and prospects. 



I) IRF Models : Basics 



The IRF models were introduce by Baxter ||2^ as Statistical Mechanical models where the 



variables are defined on the vertices of the square lattice while the interaction are defined on the 
faces. These kind of models should be distinguished from the vertex ones, where the variables 
are located on the edges and the interaction is defined on the vertices where four edges meet ||26|| . 
In certain cases there will be a deep relationship between these two types of models, given by a 
vertex-IRF map. The most interesting class of IRF models are the so called graph-IRF models 
( for a review see |^). The heights of these models are labelled by the vertices of a graph Q. 
The allowed configurations are restricted by the constraint that the lattice variables that are 
nearest neigbour in the lattice are also nearest neighbour in the graph Q. A characterization of 
Q is given by its incidence matrix A^,;, which is 1 ( reps. 0) if the heights a and b are connected 
(resp. disconnected) by a link of Q. We assume that there is at most one link connecting any 
two points. The graphs which we shall study in this paper are bipartite, which means that they 
can be partitioned into two subgraphs, say even and odd, so that any point of one subgraph is 
connected only to points of the other subgraph. A pair of variables (a, b) is said to be admissible 
if Aafi = 1. In this terminology, all the heights connected by a link of the square lattice must be 
admissible. An important example of IRF models are the RSOS models for which ^ is a ADE 
Coxeter diagram, which will be studied in detail in section IV. The Ar graph consist of r points 
labelled by a = 1, 2, . . . , r. 

In the Hamiltonian or transfer matrix formulation of IRF models, a state of the Hilbert space 
is described by the ket, 

|a >= \ao, fli, . . . , Oat > (2) 

where (oj, cti+i) is an admissible pair. There is a geometrical interpretation of the IRF states as 
paths on a Brateli diagram, which is constructed by folding the graph as in figure 1, and repeating 
the pattern along the "x-axis" |2^, A path ^ on the Bratelli diagram is a succession of points 
{^(z)}^Q such that the couple + 1)) coincides with a link of the diagram (see figure 2). 

Another important concept is that of a plaquette on the Bratelli diagram. A plaquette is the 
four-tuple, 

( d \ 

{a,b,c,d)= a ' (a, 6), (6, c), (a, d), (li, c) : admissible (3) 

V b J 

and can be identified with the elementary squares or plaquettes of the Bratelli diagram. 
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The IRF models can be defined on periodic or open chains. In this paper we shall concentrate 
on the later case. To define the dynamics of the IRF model we shall introduce the plaquette 
operator Xi, which gives the infinitesimal evolution of an IRF state in the neighbourhood of the 
z*^-site of the chain, 



a' 



Xj| . . . , aj_i, Oj, Oj+i, . . . >— ^ i? aj_i Oj+i I . . . , Oj-i, a^, Oj+i, . . . > (4) 

< \ ai J 

R{ai-i, tti, Oj+i, a^) denotes the local Hamiltonian associated to the plaquette (aj_i, a,, aj+i, a'^). If 
R{ai-i, tti, tti+i, a'i) is replaced by a Boltzmann weight then the operators are those introduced 
by Baxter in the study of integrable IRF models. In our case we are working with infinitesimal 
versions of these Bolztmann weights and, on the other hand we do not need to impose any kind 
of integrability condition, although it could be interesting to analize a possible interplay between 
the RG and integrability. 

The Hamiltonian acting on an open chain is defined as, 

N~l 

H=Y.X^ (5) 

i=l 

The time evolution produced by @ preserves the boundary heights a = and h = a^. We 
shall call 7i^^ the Hilbert space expanded by IRF states with these boundary conditions, 

'^a,b = {ko,ai, ■ • ■ ,a7v > \ \ao = a,aN = b} (6) 

Below we shall consider a generalization of this type of Hilbert spaces characterized by fixed 
boundary conditions at the ends. 

To finish this section we shall review other applications of IRF ideas in the context of Particle 
Physics, which will help us to introduce new concepts in the next section. 

It is well known the connection between integrable statistical models and factorizable S- 



matrix theories |3^. For example, the Boltzmann weights of the 6 vertex model can be conve- 
niently identified with the scattering S-matrix of the solitons of the sine-Gordon theory. This 
kind of interpretation is also possible for IRF-Bolztmann weights, which may describe the S- 
matrix of solitons (or kinks) which connect different vacua. A soliton say Sa,b{0) with rapidity 6, 
is a field configuration which connects the vacuum a at x = — oo with the vacuum b at x = +oo. 
Two solitons say Sa^bi^i) and S'f,_c(^2) can meet at the common vacuum b and after a certain time 
the middle vacuum b can turn into a new vacuum, say d. The corresponding S matrix for the 
process Sa,b{Oi)Sb,c{02) Sa,d{02)Sd,c{0i) is described by the Boltzmann weight associated to 
the plaquette (||) and rapidity 9i — 62- In this terminology the IRF state (|) can be interpreted 
as a collective state formed by N solitons connecting the vacuum and through a series of 
interpolating vacua ai, . . . , otv-i- 

There is yet another interpretation of the IRF states. If we view the graph Q as the target 
space of a discretized string, then (Q) becomes the state of an open string Sa,b with fixed boundary 
conditions at the ends. As in the case of solitons, the strings may join and split in various 
ways according to graph rules. In the rest of the paper we shall consider as equivalent the 
interpretations of the IRF states as paths on a Bratelli diagram, kinks of a field theory and 
discrete strings (see figure 3), 
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Path = Kink = String 



(7) 



II) RG of IRF-Hamiltonians: Generalities 

The basic problem we want to address is the construction of the ground state and excited states 
of the IRF Hamiltonian for large values of N. The RG method gives an approximate answer to 
this difficult problem. Wilson's strategy for the Kondo problem is to start out from small chains 
and grow them by adding site by site, while keeping a fixed, and usually large number of states, 
say m, as "representatives" of the whole chain. This method is known as the onion scheme, to 
be distinguished from the Wilson-Kadanoff blocking scheme which consists in the partition of 
the chain, or more generally the lattice, into blocks which are afterwards renormalized getting 
smaller lattices. 

The Wilsonian growth process applied to an IRF state is depicted in figure 4. A string (kink) 
with states * and a and its ends, "absorbs" a particle (vacuum) in a state h becoming a new 
string (kink) with BC's (*,&). Of course the pair (a, 6) should be admissible for the absortion 
process to be possible. The state * at the l.h.s. of the string will be kept fixed in the construction 
of "longer" strings. The strings will grow from their r.h.s. 



String Hilbert Spaces 

Let us call or more simply , the Hilbert space of the strings S which have BC's *, a at 
their ends. An example of Ti^ is given by the Hilbert space (^, where N measures the length 
of the string S. The RG method will lead to the construction of Hilbert spaces Tif which are 
imbedded into Ti^ for some N. 

The Hilbert space of the string S with an added point • on its r.h.s. will be denoted by T^f* 
and it is given by, 

"Wf* = ®h\K,b=i^h (8) 

which implies, 

dimUl' = Y.K^,diuml (9) 

b 

Proceeding as above one can construct longer strings as for example S"**. All what is needed 
is the "fusion" matrix A. For later convenience we give below the basis of the most common 
string spaces, 

nl'= {|ea®6> |A,,5 = 1} 

nl"= {\U®b®c>\Aa,t = Ab,c = l} ^^"^ 

n 

A generic Hilbert space of the form 7i^* "* will be denoted by 'Ha'^* ■ The complete Hilbert 
space of a string S plus n points •, consists in the direct sum, 
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n^^^' = ©a^f'"' (11) 

The total dimension of TC^'"'' can be computed from eq.@, 

m5,n. = dim?^^'"' = J2 K,,a,---K„.,,a^ma„ (12) 

ao,...,a„ 

where = dimTif . It is important to reahze that the sum over heights in ([T^ ) may not contain 
all the heights of the graph. For example for a bipartite graph only even or odd heights will 
appear at the right end of the string ( in this sense the strings associated to a bipartite lattice 
can be classified as even or odd). Thus if the string S is even (odd) the string plus one point S» 
will be odd (even), and S • • will again be even (odd), etc. 



String Operators 

We shall call string operators those operators O"^'"* which acting on the Hilbert space Ti"^'"*, do 
not change the height located at right hand end of the combined system S» ■ ■ ■ •. Their action 
on the basis (0) is given as follows, 

(13) 

The matrix elements of the operators O^'""* appearing in ([T3|) will be denoted by. 



<^',®h\0^'\ia®h>= of: 



*a) 



< C®d (8)010^" \^a®b®c>= oil 



(14) 



* c 
V a b 



and can be depicted as 2(n + 1)— gons, with a special vertex * from which emanate two thick 
lines representing string states labelled by C, (see figure 5). The remaining 2n thin lines connect 
admissible pairs of heights. The most important examples of string operators are given by the 
Hamiltonians H^'^* . However not all the string Hamiltonians are independent. Actually, given 
H^* and the "Boltzmann weight" R (^) one can build up the remaining Hamiltonians H^'"'* for 
n > 2. The first member of the later family , namely , has to be given independently, but 
quite paradoxically it plays little role in the construction. 

As as example we give below the matrix representation of H^", 
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h: 



/ e d 

* c 
\ a b 



(15) 



i/g I * 6 1 6b,d Ab,c + Sa,e ^5.,^ ^ 



a c 
b 



This eq. is depicted in figure 6, where we show also the construction of H^'^'. 



The RG-operation 

The key point of the RG method is the construction of the RG-operator T that truncates the 
Hilbert space Ti.^* into Ti.^ , where S' represents a string with one more site than the string S, 
i.e. 



T : 



n 



s' 



The matrix representation of T and its hermitean conjugate T"^ are given by, 



(16) 



(17) 



where 



According to ([l6|) T is a m^/ x ms» matrix. Except for the first RG-operations we shall 
always keep the same number of states describing the renormalized system, i.e. m = ms = ms'- 
Both T and can be depicted as triangles with the special vertex *, which is the origin of 
two thick edges which symbolize the old and new (renormalized) strings ( see figure 7). The 
truncation operator must satisfy the equation. 



TTt = 1 



(19) 



which guarantees that T"'"T is a projection operator which maps Ti^* into a subspace which is 
isomorphic to Ti^ . 

Eq. ([l9|) reads in components (see figure 8), 



a I 'ib \ * b j ^'''^b 



(20) 



Given the operators T and we can renormalized any operator O^'"'* down to an operator 
qs ,(«-!)• -j^y jY^gans of the equation. 
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QS',in-l). ^ J. (^S,n. yt^ > 1) (21) 

which in mathematical terms expresses the commutativity of the following diagram, 

^5',(n-l). O^^^D- ^5',(n-l). 

Tt i ] T (22) 



In eqs.(^) and ( |2^ ) the operators T and T'l" act trivially on the points beyond the closest 
one to the string S. As an example we give below the renormalization of O^'* and O"^'^* (see 
figures 9 and 10 ), 



* "1^' .^ (23) 



( 

\ a 



d a 

In summary we have presented in this section a formalism to deal with the renormalization of 
generic IRF Hamiltonians. In the next section we shall explain the DMRG algorithm to construct 
the truncation operator T, which will then allow us to carry out explicit computations. 



Ill) The IRF-DMRG algorithm 



The "standard" RG method to construct the operator T applied to IRF models consists in the 
following two steps: i) diagonalization of the Hamiltonian Ti,^* and ii) projection to its lowest 
energy states. This algorithm treats the system S» as isolated from the rest of points which 
one adds in posteriori RG steps. In other words, the height associated to the point • in S» is 
fixed to a given value. Imposing fixed boundary conditions at the ends of the blocks in the RG 
method always leads to bad results. Instead one should consider a combination of B.C.'s as in 



3^ , or impose open B.C.'s as in [Q. 



The DMRG is a way to take care of the influence or correlations of those points that have 
not yet been added to the block. There are various DMRG algorithms: infinite system method, 
finite system method, etc. We shall give in this paper the IRF version of the infinite system 
algorithm, which is based on the superblock formed by a string S, a point • and another string 
5''^, which is the mirror image or reflection of the string 5* ( see figure 11). The dynamics of 
the "super-string" S • involves all allowed heights at the middle point •, and in that way 
one is not commited to a particular B.C. on •. There is an appealing electrostatic analogy to 
understand the role of S^. Let us recall the mirror image method which is used to impose 
Neumann (open) B.C.'s on the electrostatic potential. In this sense the mirror string seems 
to play a similar role, i.e. that of imposing open B.C.'s on •. 

A basis of the Hilbert space of the super-string 5* • 5*^ is given by. 



7 



n''' ={|e,®6®77,> II A,,, = A,,, = 1} (25) 

The Hamiltonian which generates the dynamics of the states belonging to Ti^*^"^ can be 
obtained using the methods of the last section and it reads (recall figure 6), 



a' b' d ^ 




{ ' 




* * 


* b 




a b c J 








( b' \ 






( c 



(26) 



b I V 



* b 
c 



Now we diagonalize this Hamiltonian and select its ground state which is called the target 
state and can be written as, 

1^0 >= EE^C,r,(a,6,c) |e®6®r/> (27) 

a,b,c 5,77 

The mirror string plays an auxiliary role in the construction and we should get rid of 
it. The DMRG proposal is to construct the reduced density matrix p^* of the subsystem S» by 
tracing over the states in S^, 



Tr^sfl |?/'o >< ^ol 



(28) 



In the above trace we shall set the height of the middle point the same for both the ket and the 
bras, so that the matrix representation of p^* will be given by. 




= im^€>^(«'^ c) ^l,^^{a,b, c) 

c -q 

A normalized ground state ( pTj) yields a properly normalized density matrix. 



(29) 



a,b 5 




(30) 



The next step is to diagonalize the matrix (^) in the Hilbert space Tif*, for every value of 
b, keeping the first m eigenstates with highest eigenvalue. These states are the most probable 
ones to contribute to the ground state of the super-string. Finally the matrix T ([16|) is given by 
these m column vectors. 

Eq.(p9|) is very similar to Baxter's definition of the corner transfer matrix (CTM) for IRF 
models 0], in the sense that one traces over the degrees of freedom of half of the system while 
keeping the height located at the edge of the "cut" fixed. This relation between the DMRG and 



the CTM has already been pointed out in |jTT[, and we expect it to hold also for the IRF-DMRG. 
This ends our presentation of the IRF-DMRG. 
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IV) The IRF-DMRG at work 



We shall apply below the formalism developed in the last two sections to study IRF models 
that can be obtained by means of vertex-IRF maps of vertex Hamiltonians. This map will be 
explained in detail in section V for the case of the SOS models. 



SOS model (S=l/2) 

The spin chain Hamiltonian of the Heisenberg model with spin 1/2 reads, 



^=^E(^^^m + l) (31) 



where ai are Pauli matrices acting at the i^^ site of the chain. The choice of (^) is motivated 
by the fact that \ (ajcrj+i + 1) is the permutation operator acting at the sites i and z + 1. The 



model defined by ( plD is equivalent to an IRF model whose graph, denoted by Aqo, consists in 
the semi-infinite chain of fig. 12. The heights j = 0, 1/2, 1, . . ., that label the points of the graph 
Aoo, are in one-to-one correspondence with the irreps of the group SU{2). According to fig. 12 
the incidence matrix of satisfies 



A,,y = l ^ \j-f\ = 1/2 (32) 

The IRF-Hilbert space associated to a chain with N-sites is given by the direct sum ( recall 
(0)), 

= Q^nf (33) 

where TC^ is the IRF-Hilbert space of all the states with total spin j, 



nf = {\jo,ju...,jN> \\ 30 = 0, jN = j, \j^-j^+l\ = 1/2 for i = 0, . . . , iV - 1 } (34) 



According to ([3^ ) the height -k should be identified with the identity irrep (i.e. * = 0). Since 
the graph A^o contains an infinite number of heights this model is an unrestricted IRF model 
called solid-on-solid model (SOS). This implies in particular that the Bratelli diagram consists of 
a pyramid of infinite height as one moves from the origin of the diagram (i.e. * point = vacuum 
representation) to the right hand side ( fig. 13). 

The dimension of the Hilbert space (0) is given by, 

/N 



This formula can be compared with the number of vertex states of the standard formulation 
of the Heisenberg model, with a fixed value of the third component of the spin ( see section 
V), 

dimViI= f 1 (36) 
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For N even (odd) the ground state will belong to the Hilbert space with j = (1/2). As can 
be seen from eqs.(p5D and (^) it is more efficient, for numerical purposes, to look for the ground 
state of the Heisenberg model in the IRF subspaces than in the vertex ones ( see table 1). 



N 


dim7i:jlo(IRF) 


dim?-^ft^Q (vertex) 


4 


2 


6 


10 


42 


252 


20 


16 796 


184 756 


24 


208 012 


2 704 156 


N » 1 


~ 1. 59562^ /Ar3/2 


~ 0.79782^ /iVi/2 



Table 1 



From (^5]) and (|36| ) we get the relation, 

dimViY_ 



dimTij^o 



(37) 



which is a numerical version of the eq.(|T]) and shows that the difference between the vertex and 
IRF formulations persists in the thermodynamic limit. 

The constraints (|3^) imply that there are only 6 different "Boltzmann" weights whose values 
are given in table 2. 



d 




( \ 


a c 


R 


a c 


b 




[ b J 


J 

J ± 1/2 J T 1/2 




1 


J 






J ± 1/2 
J 3 
3 ± 1/2 


2j+l 


3 ± 1/2 
3 3 
3 T 1/2 


V2i(2j+2) 
2j+l 



Table 2: The SOS (S=l/2) Hamihonian. 



Numerical Results 



In table 3 we present the data for the ground state energy of the Heisenberg Hamiltonian (pi]) 
for chains of length = 6 up to 24. In this IRF-DMRG computation we keep m=12 states, 
which is a rather modest number of states, while for the vertex-DMRG we keep a maximum of 
12 states , since it is impossible to fix the number of m for vertex-DMRG, due to the degeneracy 
based on SU(2) symmetry. This is one of the advantages of the IRF-DMRG as compared with 
the vertex-DMRG. It is clear that at equal number of retained states, the IRF method should 
give better results than the vertex method. This expectation is conffimed in table 3. 
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N 


Exact 


IRF-DMRG 


Vertex DMRG 


6 


-2.4871542677758 


-2.4871542677758 * 


-2.4871542677758 * 


8 


-3.2498651973757 


-3.2498651973757 * 


-3.2498651973757 * 


10 


-4.0160704145657 


-4.0160704145657 * 


-4.0160641768009 


12 


-4.7841812656810 


-4.7841812656810 * 


-4.7835746471807 


14 


-5.5534493237243 


-5.5534493236562 


-5.5527041895949 


16 


-6.3234742911502 


-6.3234742887647 


-6.3246392674964 


18 


-7.0940221370730 


-7.0940221268443 


-7.0953356164320 


20 


-7.8649466687979 


-7.8649466378157 


-7.8663701424885 


22 


-8.6361517519671 


-8.6361516790424 


-8.6375643759289 


24 


-9.4075715208191 


-9.4075713713902 


-9.4089821742785 



Table 3: Ground state energy of the Hamiltonian (|3TD. The data followed by "*" are exact. 

If we increase the number m of states retained, the results converge exponentially fast both 
in the vertex-DMRG and IRF-DMRG methods ( fig. 14). 



RSOS models 

An interesting generalization of the spin 1/2 Heisenberg chain is provided by the XXZ Hamil- 



tonian with boundary terms 



H 



xxz 



W-1 



El ^^^^ I ^Y^Y I ^ ^ ^ 



q-q 



-1 



'i ^i+i 



Nj 



1=1 



(38) 



where q = e^"* is a phase. This Hamiltonian has very interesting properties: 



The eigenenergies of the N=2M site XXZ chain (P3) coincide with those a M-site self-dual 
Q-state Potts model with Q = q + q~^ = 2 cos'j p7| . 

Invariance under the action of the quantum group SU{2)q [pH] . 



Using q- group theory one can map the vertex Hamiltonian ([38| ) into a RSOS Hamiltonian 
whose graph is given by the Coxeter diagram ( see figure 15) |29(| . 



jjxxz -g Qj-j^icg^i gj^(^ ^Qj- J = _E_ belongs to the universality class of the minimal 
CFT's pOl with a value of the central charge given by, 
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r(r + 1) 



(39) 



We shall study below the RSOS version of the vertex Hamiltonian (|38|). A way to arrive to this 
version consists in writing (|38|) as follows 



Hxxz ^ j^N-i _ (40) 

where are the Temperley-Lieb- Jones (TLJ) operators which act at the positions i^^ and {i+iy^ 
of the chain and satisfy the TLJ algebra 
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= (g + g ^)ei 

Ci (ij ~ (ij (^ij \i J I ^ 2 

In the vertex basis the TLJ operator Cj can be written as follows, 



(41) 



ei = li 



h-i ® 



(S) 1 



i+2 



(42) 



^ \ 

q-^ -1 
-1 g 

V y 

The existence of a vertex-IRF map, and the fact that the operators Cj commute with the 
action of SU{2)g imply that they can be given a representation on the RSOS-Hilbert spaces of 
the face model defined by the graph A^, 



6j I . . . , Clj, 



...>=E 



e I aj_i Oj+i j I . . . , tti-i, a'j, Oj+i, . . . > 



d 



e \ a c = 



tn 



(43) 



(44) 



where are the components of the Perron-Frobenius vector of the incidence matrix of the graph 
Ar, which are given by 



tn = sin 



na 



r + 1 



, a = l,2,...,r 



Notice that ta satisfies, 



ta+l + ta-l = 2C0S 



TT 



r + 1 



(45) 



(46) 



Recalling that the incidence matrix satisfies in this case 



A. 



1 



|a-6| = 1 



(47) 



one gets that there are only 6 types of "Boltzmann weights", whose expression, given in table 
4, can be computed using eqs (|40|) and 
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d 




/ 1 \ 

( d \ 


a c 


R 


a c 


b 




{ b 1 


a 






a ± 1 a =F 1 




COS7 
2 


a 






a±l 






a a 
a±l 


COSf tail 
2 ta 


a±l 
a a 
a =F 1 




\/ta + lta-l 
ta 



Table 4: The RSOS Hamiltonian for the face model A^. 

In the limit where r ^ oo the RSOS model becomes equivalent to the SOS model 
with graph A^ studied previously. One can check that the R-matrix given in table 4 is, up to a 
constant and a change of basis, the same as the R-matrix given in table 2, with the identification 
a = 2j + 1. 



Numerical Results 

In table 5 we give the ground state energy per 2 sites Eq{M)/M{M = N/2) of the XXZ model 
(PSp, which coincides with the ground state energy per site of the corresponding Potts model. 
This table should be compared with table lb in reference |3^, which was obtained using the 
Bethe ansatz. The authors of give their results up to 6 decimals ( ours is 9) and the agreement 
in the energies holds until the 6*^ digit. The number of states retained in our computation is 
m = 160. 

Using the IRF-DMRG data we have computed the finite size corrections to the ground state 
energy, which are governed by the formula fi^, 



Eo{M)/M = eoo + § - ^ + o(M-^) (48) 

where Coo and f^o are, respectively, the bulk and surface energy per site. <^ can be identified with 
the spin wave velocity and it is given for the Potts model by, 

c = ^ (49) 

27 

We have used Sach's formula to get the values of the central charge c |^ . 

The outcome of this computation is that the IRF-DMRG method reproduces rather accu- 
rately the results obtained using the Bethe ansatz. This supports the hypothesis that the DMRG 
is in fact an exact numerical RG method. 
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N/2 


r=3 (Q=2) 


r=5 (Q=3) 


r=7 (Q=3.414) 


r=inf (Q=4) 


2 


-1.320 899 500 


-1.478 675 250 


-1.537 532 848 


-1.616 025 403 


4 


-1.459 958 153 


-1.580 754 400 


-1.626 304 313 


-1.687 466 299 


8 


-1.532 472 880 


-1.636 075 843 


-1.675 237 896 


-1.727 934 286 


16 


-1.569 617 420 


-1.665 066 291 


-1.701 130 752 


-1.749 664 452 


32 


-1.588 433 915 


-1.679 941 590 


-1.714 488 941 


-1.760 960 537 


64 


-1.597 906 426 


-1.687 482 175 


-1.721 280 643 


-1.766 733 433 


128 


-1.602 659 177 


-1.691 279 439 


-1.724 706 257 


-1.769 649 934 


256 


-1.605 039 733 


-1.693 185 000 


-1.726 426 767 


-1.771 116 470 


512 


-1.606 231 063 


-1.694 139 540 


-1.727 288 989 


-1.771 851 868 




-1.607 423 097 


-1.695 095 264 


-1.728 152 544 


-1.772 588 719 


(AB3Q) 


-1.607 423 


-1.695 095 


-1.728 152 


-1.772 588 


foo 


0.610 501 838 


0.489 636 433 


0.442 486 891 


0.377 747 124 


(AB3Q) 


0.610 502 


0.489 637 


0.442 487 


0.377 649 


c 


0.499 942 


0.798 817 


0.893 150 


1.038 18 


(AB3Q) 


0.500 00(1) 


0.799 9(2) 


0.89(3) 


0.99(2) 


(exact) 


0.5 


0.8 


0.892 857 


1 



Table 5: Ground state energy per 2-sites of the RSOS chain. We also give the results of ref ||2TI| 

(AB^Q) for Coo, foo and c. 



V) Vertex-IRF Map 

We shall illustrate this map in the case of a spin-s chain whose dynamics is dictated by a 
rotational invariant Hamiltonian. The Hilbert space of the spin-s chain with N sites will be 
denoted by V^''^ and consists in the tensor product of N copies of the vector space Vg = C^'^"''^, 
where acts the local spin operators S, ( i = 1, . . . , A^). The vertex-IRF map is based on the 
tensor product decomposition of the space V^'** into its irreducible components. 



Vertex Hilbert Spaces IRF Hilbert Spaces 

Using the Clebsch-Gordan decomposition of tensor product of irreps of SU (2) one can write, 

yN,s^ ^ nf^'^Vj (50) 
0<i<2sAf 

where Tif'^ is the generalization of the IRF Hilbert space (|3^) to the spin s case. The heights 
tti G that label the IRF states are subject to the following constraints, 

ao = 
01 = 5" 

\S — ai\ < Oj+i < 15* + Oil, i = l,...,N — l 
aN = j 

The dimension of TCf'^ is given by the number of times the spin-j irrep appears in the CG- 

N 

decomposition of the tensor product ■ ■ ■ ®s. 
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dimHf' = multiplicity of in V^'" (52) 
diniTi^''* can be computed using the following formula, 

dimTi^'"* = dimVj^'** — dimVj^i (53) 

where VjY'* denotes the subspace of V^'"* with a fixed value of the third component of the spin. 
Eq. ( p^D says that the highest weights with total spin j are given by the states with spin = j 
minus the ones that can be obtained from Vjl+i by the lowering operator S~ . The counting 
of states with a fixed value of is easily done using the "Bethe method" of starting with the 
ferromagnetic state with all the spins up and lowering the spin. Below we give the formulae for 
s=l/2 and 1. 



N,s=l/2 



dimV, 



dimVi!'^=^ = eE"^^^^" 




(54) 



s"" + k 



where the symbol [x] appearing in the upper limit of the sum denotes the integer part of x. 



The relation between the vertex basis of the spaces V^''*, Vj and the IRF basis of Tif'^ can be 
obtained using the Clebsch-Gordan coefficients as follows, 



.,mjv 



s ai 
mi ni 



a,N-2 S CtAf-l 
O'N-2 nj^^i 



where rij = mi + ■ ■ ■ + m^ 



m 



riN 



e(a) 



ai s 02 
rii 1712 n2 

a^^i s 
riN-i rriN 



a2 s 



as 

n3 



(55) 



N 



rrii 



and a denote the IRF labels which satisfy 



conditions (|51|). A graphical representation of eq.(|55D is given in Fig. 16. The at the upper 
left of the diagram can be identified with the * symbol introduced in section II. The vertex-IRF 
map, as defined by eq.(^), is nothing but a change of basis from vertex variables to IRF ones 
which achieves the factorization of the SU (2) symmetry. 



Vertex-Hamiltonians IRF-Hamiltonians 

The most general form of a rotational invariant Hamiltonian H acting in V^'*, i.e. 



N 



[H,S] = 0, S = J2S., (56) 
1=1 

which is translational invariant and contains only nearest neaghbours couplings is given by, 

N-l 2s 



(57) 



i=l r=0 
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where is a set of 2s+l coupling constants (oq can be put equal to zero since it multiplies the 
identity operator). Since (|57| ) commutes with SU{2) it means that its action affects only the 
IRF spaces Hf. Using group theoretical methods we get, 



Hasi) = j:l~,' r( )e(---,a.-i,a>m,---) (58) 

where the IRF "weights" R can be computed in terms of the coupling constants and the 
6j-symbols as follows (the details of this computation will be given elsewhere), 

Ja,-i < \= I ' ' M I ^ ' ^ I (59) 



Af = E'loarX,^ (60) 
Xj = |j(j + 1) - s{s + 1) 

As an example we may choose H to be the sum of all the permutation operators between 
nearest neighbours, in which case Aj^ turns out to be a sign factor, 

H = Y: Pi,i+i =^ = (-1)'^^^' (61) 

i 

The IRF hamiltonian corresponding to (|6TD is given by, 

( j _ ^_^yi-i+ai+i-ai-a[ I ^ '^i-l '^i 1 ^g2) 

\ flj J \ s ttj+i a'- J 

The s=l/2 Heisenberg Hamiltonian ( pTj ) is precisely of the form ( |6lD so that table 2 can 
be derived from (|62|). In a subsequent publication we shall use eq.(^9|) to study higher spin 
Heisenberg chains in the IRF formalism. 



Vertex-DMRG ^ IRF-DMRG 

The DMRG algorithm to renormalize a block B plus one point into a new block B', is based 
on the superblock B QB^, where B^ is the reflection of the block B ( we use to distinguish 
vertex-points from IRF-points which were denoted above by • ). The main steps of the vertex- 
DMRG are: 

• Diagonalization of the superblock Hamiltonian to find the ground state {ip >. 

• Construction of the reduced density matrix by tracing over the states in QB^, 

p^® = TrQBi^ \iJo><H (63) 

• Diagonalization of p^® to find the eigenvalues Wa and eigenvectors |m" >. Discard all but 
the largest m eigenvalues and associated eigenvectors. 



16 



• Construct the operator T using the truncated eigenvectors >. 

• Renormahze all the operators using the analog of eq. (^Tj) in the vertex case. 

• Repeat the process for the new block B'. 

These set of rules define the vertex-DMRG algorithm, which applies directly to systems where 
the lattice variables associated to the points © are not subject to constraints except perhaps for 
conservation laws like total spin, charge, etc. Most of the Hamiltonians in Condensed Matter 
or Stat. Mech. are of this form. If the vertex Hamiltonian happens to have a continuous 
symmetry, then the factorization of that symmetry would lead naturally to an IRF model, 
whose renormalization can be studied using the IRF-DMRG method. 

The relation between the vertex-DMRG and the IRF-DMRG algorithms, presented in section 
III, is illustrated diagrammatically in fig. 17: the height a (resp. b) labels the different irreps 
(f.ex. the total spins of the spin chains studied above) that appear in the tensor product of 
all the irreps contained in the block B ( resp. B^). The intermediate height j is obtained 
tensor ing a (resp. b) with the irrep carried by the vertex 0, which in the case of the spin chains 
is a spin-s irrep. Finally, we must tensor j ® j and pick up the identity irrep. Let us now find 
the analytic relation between the vertex and IRF density matrices. We shall call ( resp. V^'^ 
) the Hilbert space associated to the block B ( resp. B^). The tensor product decomposition 
(^) becomes in this case 



V^ = Ea^f®V; (64) 
Using this decomposition the Hilbert space of the superblock B QB^ becomes, 

yBQQB^ = ^ T^f K ® K ® K ® H ® nf (65) 

a,b 

The ground state of any rotational invariant Hamiltonian acting in this superblock can be 
written in the basis of ([65| ) as follows (see fig. 18), 



1^0 >= E l^a ® ma® mi® 1712® nib® rib > i)^^,riM,j, b) 



(66) 



j J 

ma + mi mb + m2 



a s ] 

ma mi ma + mi 



s b j 

7712 mb m2 + mb 



where i'^a,vbi^jjy ^) I-^F wave function of the ground state >. 

The density matrix p^® can be obtained from (|63|). Using the properties of the CG coefficients 
we get. 



P 



B& 



m„ 



mi>< i'^, 



m„ 



(67) 



a' 1 n a 

2j + l y/ia 



a s J 

uia mi rUa + mi 



m„ 



s 3 
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where 



a 

* 3 
a 



(68) 



coincides with the definition of the IRF-DM given in (^). The relation (|67D between the vertex 
and IRF density matrices can be finally written as, 



P 



B& 



E 



r2j + 



p^ 



(69) 



The factor l/(2j + 1) takes care of the degeneracy of the irrep Vj in the CG decomposition 
a® s ^ i, and guarantees the correct normalization conditions of both density matrices. 



VII) Conclusions and Perspectives 

We have generalized in this paper the DMRG method to Id Hamiltonians of IRF type and 
showed, in the examples of the spin 1/2 SOS and RSOS models, that it gives very accurate 
results. Our method is equivalent, by means of a vertex-IRF map, to the standard DMRG 
method formulated by White for Hamiltonians of the vertex type. This map consists in the 
factorization of the symmetry group of the vertex theory. This factorization has numerical and 
conceptual advantages. From a numerical point of view one needs to keep a smaller number of 
states in the IRF-DMRG in order to achieve the same accuracy as in the vertex-DMRG. The 
degeneracy of the eigenvalues of the vertex formulation, due to the symmetry, is absent in the 
IRF case, which makes the numerical analysis more compact and stable. Conceptually the IRF- 
DMRG is also very appealing since it employs tools and techniques well known in Statistical 
Mechanics, Integrable Systems, Multi-Matrix Algebras and Conformal Field Theory. Thus the 
IRF states can be seen as paths of a Bratelli diagram, kinks of a Theory of Solitons, discretized 
strings and conformal blocks in CFT. The formalism we have developed allow us to apply the 
DMRG method to IRF states in a very natural way. 

Let us mention some of the lines of research which we believe deserve further study. 



Higher Spin and Ferromagnetic Spin Chains: In section V we have presented the 
necessary tools to study higher spin chains. A particular interesting case is the spin 1 
chain , which has a rich phase diagram. We shall show in a subsequent publication that 
the string order parameter of den Nijs and Rommelse |4^, which is used to characterize 
the Haldane phase, adopts a particular simple form when written in IRF variables. In fact 
the IRF states constitute a complete and orthonormal basis of valence bond states. In 
particular the AKLT state [Q, which is a pure Haldane state, is simply a straight path in 
the Bratelli diagram of the spin 1 Heisenberg chain. 

The IRF-DMRG is also very promissing for the study of ferromagnetic systems, which 
seems to display a rich phase structure in the presence of magnetic fields . The vertex- 
DMRG method applied to ferromagnetic systems encounters the difficulty that the ground 
state has a huge degeneracy. As we have shown in this paper, the IRF-DMRG eliminates 
this degeneracy, avoiding the complications arisen from that fact. 
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t-J and Hubbard models: The vertex-IRF map can be straighforwardly applied to 
these models yielding an IRF formulation where the spins form valence bonds. The IRF 
heights are now given by the couple (spin (j) , charge (q)). In the case of Hubbard model 
the symmetry group is given by 5*0(4) and it contains, in addition to the rotational group, 
the group of pseudo spin rotations. The factorization of this larger group should reduce 
considerably the dimension of the Hilbert Spaces. 

Ladders: These systems, which have received considerably attention in the last 2 years, 
consist of a finite number of coupled chains, with very interesting properties ( for a review 
see [0). For spin ladders with a few number of chains it is rather simple to obtain their 
IRF version simply by taking the tensor product of the irreps located on the rungs and 
performing afterwards their tensor product along the chains. This procedure imitates the 



strong coupling analysis applied to these kind of systems The IRF models so obtained 
has more than one link connecting different heights, and so one has to generalize slightly 
the construction of this paper. A similar multiplicity phenomena occurs in the theory of 
solitons [0. The IRF formulation of ladders could be useful to clarify the relationship 
between their phases and those appearing in higher spin chains. 

• Higher Dimensions The DMRG philosophy is not confined to Id, but the standard 
DMRG algorithms proposed so far are one dimensional, despite of some 2d applications to 
finite clusters ||10[. The IRF-DMRG strengthen this point of view, since in particular the 
vertex-IRF map is a one-dimensional operation. We should perhaps say that the vertex- 
IRF map is really adimensional because the tensor product operation does not impose 
any particular geometry or dimension. In connection with this problem it may be useful 
to realize that the vertex-IRF map is a duality transformation similar to the Krammers- 
Wannier duality or the Jordan- Wigner transformation. This interpretation may serve as a 
guide to construct higher dimensional vertex-IRF maps. 

There are still many more topics to be considered in connection with the DMRG. The DMRG 
method has arisen as a numerical tool specially well adapted to Id systems, but in our opinion 
its importance goes beyond its numerical success. There are still some fundamental questions 
whose solution we would like to know. It is perhaps not exagerate to say that new and radical 
developments connected with the DMRG are likely to happen in the near future. 
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Figure Captions 



Fig.l.- Bratelli diagram associated to the Coxeter graph A^. 

Fig. 2.- In dark it is shown a path on the Bratelh diagram of Fig.l. 

Fig. 3.- The string S'^.^^ as a representative of the class of all paths on a Bratelli diagram starting 
at * and ending at a. 

Fig. 4.- A string S^:,a absorbs a point • which carries an allowed state b, becoming a new string 
Fig. 5.- Diagrammatic representation of the string operators (p!4D. 

Fig. 6.- Top: diagrammatic representation of the equation (plsl) . Bottom: Diagrammatic recon- 
struction of the Hamiltonian H^'^'. 

Fig. 7.- Diagrams of the T and T"!" operators ([TTD ( T for truncation and for triangle). 

Fig. 8.- The normalization condition (^) interpreted as a kind of annihilation process of triangles. 

Fig. 9.- Eg . (p3D in pictorical form. 

Fig. 10.- The renormalization of O^'"^* ( see eq.(^)). From figs. 8. 9 and 10 we see that the RG 
procedure is a kind of sewing or gluing construction involving triangles, plaquettes and higher 
n-gons. 

Fig. 11.- The "super-string" configuration that leads to the infinite system IRF-DMRG algorithm. 
Fig. 12.- Coxeter graph A^o. 

Fig. 13.- Bratelli diagram built up using the Coxeter graph Aoo- 

Fig. 14.- Plot of the deviation of the IRF-DMRG ground state energy of a s=l/2 chain with 512 
sites, as a function of the number of states retained m. 
Fig. 15.- Coxeter garph Ar. 

Fig. 16.- Graphical representation of the vertex-IRF map. Notice that the IRF points • and the 
vertex points belong to lattices which are dual one another. Indeed the vertex-IRF map is a 
kind of duality transformation. 

Fig. 17.- The vertex-IRF map that relates the vertex-DMRG and the IRF-DMRG algorithms. 
Fig. 18.- Here we show the CG decompositions involved in fig. 17. 
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